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I. INTRODUCTION 


An important physical phenomenon which has many practi- 
cal consequences in technology, ranging from fuel injection 
sprays to ink jet printers, is the breakdown of a circular 
jet of fluid issuing into a quiescent environment 

This problem has been extensively studied previously 
by Rayleigh and others, [Refs. 1, 2 and 3], for steady uni- 
form jets issuing from an orifice , but’ only recently have 
the effects of jet velocity profiles or other initial condi- 
tions on the flow been investigated, [Refs. 4, 5 and 6]. 

Theoretical investigators have typically assumed jet 
breakdown to be physically similar to transition from lami- 
nar us turbulent flow and have examined the phenomena using 
classical stability theory. 

Rayleigh , [Ref. 1], studied the inviscid stability of 
Poisseuille flow and concluded that the flow was stable to 
infinitesimal disturbances. This work attracted the interest 
of many investigators and sixty years later Tatsumi [Ref. 2] 
published his work related to the study of the velocity dis- 
tribution in the axisymmetric laminar inlet flow through a 
circular tube, under the assumption of ' almost similarity ' 
of velocity profiles . He obtained fairly good agreement of 


theoretical with experimental velocity distribution curves 
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in a rather limited region near the entrance where the ap- 
proximation was valid . 

Batchelor and Gill, [Ref. 3], have analyzed the stabili- 
ty of steady axisymmetric jets considering the typical small 
disturbance as a Fourier component with sinusoidal dependen- 
ee on both axes in ecylindrical coordinates . 

Crow and Champagne, [Ref. 4] , have employed a loud- 
speaker to generate a wave train on a turbulent jet and 
found that the phase velocity of the waves could be descri- 
bed by applying linear theory of temporally growing insta- 
bility to a 'top hat' velocity profile. 

Salwen , Chester and Grosch, [Ref. 5], have studied the 
stability of Poisseuille flow in a pipe of circular eross 
section with axisymmetric disturbances . They formulated the 
problem by expanding the perturbation velocity and pressure 
in a complete set of orthonormal funetions which satisfied 
the boundary conditions . 

In a recent work by Anderson and Bejan, [Ref. 6], the 
authors presented a linear stability analysis of the large 
scale structure of a round jet surrounded by an annular she- 
ar layer . Their study was limited to the developing region 
near the jet nozzle at very large Reynolds number ( Re--o), 
They examined the radial dependence of the amplitude of gro- 
wing disturbances in order to illustrate the extent to which 


the disturbance penetrates into the jet and its surroundings 


12 


and found that the region influenced by a disturbance was 
directly proportional to the wavelength of the disturbance. 
They also pointed out that amplified disturbances exhibit a 
phase lag across the shear layer , which may account for the 
spade~like structures evident in flow visualizations of tur- 
bulent jets. 

In the present work we seek to investigate experimental- 
ly the effects of small sinusoidal perturbations on the 
breakdown of circular jets issuing from long tubes , which 
permit some control of their mean velocity profiles , for 


both fully developed laminar and turbulent flows. 
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II. EXPERIMENTAL APPARATUS AND PROCEDURE 


A. DESCRIPTION OF EXPERIMENTAL APPARATUS 

The experiment was carried out in a vertical falling 
jet produced at the outlet of a 70.0 ineh long length of 
0.1875 ineh inside diameter stainless steel tube . The inlet 
of the tube was fed from a phenum chamber machined into an 
aluminum block , which also contained a 75 Watt loudspeaker 
driver unit mounted properly to this block, used to introdu- 
ee the pressure perturbations . A schematic diagram of the 
experimental apparatus is illustrated in Fig. 1 . Details of 
the aluminum block and adapter are shown in Fig. 2. 

Pressure fluetuations as well as mean pressure within 
the phenum chamber was monitored by the Statham differential 
pressure transducer shown in Fig. 1 and a mereury manometer. 
Mean water flow was controlled with a needle valve and mea- 
sured with a 0.6 gpm rotameter as shown in Figs Lowi order 
to minimize effects of geometry at the exit end of the steel 
tube, it was very carefully machined to ensure sharp edges 
and an absence of burrs. 

The loudspeaker driver was excited with a Wavetee audio 
oscillator capable of producing sine waves , as well as 
triangular, unit funetions and ramp funetions from 1 Hz to 


10 MHz, amplified by a Hewlett Packard power amplifier. 
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Wavetec H.F VCG Generator , Model 142 
Hewlett Packard Power Amplifier , Model 467A 
Tektronix Dual Beam Oscilloscope , Type 551 
High Frequency Stroboscope , Model 839 
Cimeon Digital Multimeter DMM 51 

Hewlett Packard RMS Voltmeter , Model 3400A 
Intronics Power Supply +15 VDC 

Calibration Bridge 

110 VAC Outlet 

10 Rotameter 0.6 gpm 

11 Statham Pressure Transducer PM131 TC 

12 University Sound Driver Unit [D-75 

13 Mercury Manometer 
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Figure 1. Schematic Diagram of Experimental Apparatus 
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Figure 2. Aluminum Block and Loudspeaker Mounting 
Adapter Details 
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Excitation voltage as well as the output of the plenum cham- 
ber pressure transducer were monitored on a two channel Tek- 
tronix oscilloscope. 

Visual observations of the jet breakdown process was 
made with the aid of a high frequency LED stroboscope, deve- 
loped at the Max-Planck Institute for Fluids Research? which 
produces a constant mean light output from 100 Hz to 1 MHz. 
This device not only may be set for any phase angle but also 
contains circuitry which produces’ automatic phase angle 
sweep at an adjustable rate which permits periodic phenomena 
to be visualized as they develop. Unfortunately the light 
output was not sufficient for photographic recording. 

Photographic recording of the observed phenomena was 
made with the illumination from a General Radio Microflash 
unit, which altough eifichted than the LED stroboscope, lac- 
ked the sophisticated triggering circuitry to permit control 
of the phase angle . A 4x5 inch view camera fitted with a 
210 mm objective was used to photograph on Kodak Tri-X film 
which was force developed to achieve an effective film speed 
of ISO 3200 . 

An overall view of the experimental apparatus concern- 


ing and the photographic equipment used is shown in Fig. 3. 


IThe author extends his oraeieude to Dr. G.E.A Meier 
who made the instrument available for this work. 


Ea 





Figure 3. Photo of Overall Experiment Setup 
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B. EXPERIMENTAL PROCEDURE 

The first part of the experiment was devoted to calibra- 
tion and determination of the steady flow behavior of the 
jet without the introduction of external perturbations. 
As the flow rate was increased a laminar flow column was ob- 
served to develope which converged and began to breakdown 
into a turbulent flow at 100 to 125 diameters downstream of 
of the tube eat corresponding to Reynolds numbers of 880 to 
2200 . 

Increasing the Reynolds number above 2200 resulted in 
a turbulent jet which was sudden and easily identified with 
the naked eye . Breakdown of the turbulent jet was observed 
with the high frequency stroboscope and was observed to ran- 
ge from 595 diameters downstream at Reynolds number of 2288 
to 25 diameters as the Reynolds number was increased to 4400 

Based on these preliminary results , a subcritical Rey- 
nolds number of 1760 was adopted for laminar flows and 2640 
for turbulent flows andthe effects of small perturbations 
on the flow at these two Reynolds numbers was investigated 
over the frequency range of 0 to 2000 Hz using the high 
frequency stroboscope to observe the flow. 

Based on these observations , photographic records were 
made at selected frequencies at which significant characte- 
ristie phenomena occured . Because of the difficulty in 


accurately describing the pressure fluctuations actually 
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produced in the phenum chamber, oscilloscope traces of these 
pressure variations as well as the driving signal were also 
recorded. These results are presented and discussed in Chap- 
ter IV. A disecusston of the effeets of plenum pressure vari- 
ation on the flow at the exit of a long tube is presented in 
Chapter III .An analysis of the dynamic response of the 
pressure eneneaucer and its calibration curves may be found 


in Appendices C and D. 
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Ili. THEORETICAL ANALYSIS 


A. BREAKDOWN OF LAMINAR PIPE FLOW 

One of the most important problems in fluid mechanics, 
one which has attracted the interest of investigators for 
many years , is that of the transition from laminar’ to tur- 
bulent flow. 

The theoretical analysis of the stability of laminar 
flow is based upon the assumption that turbulence is produ- 
ced by the growth of what are initially small perturbations 
or disturbances of the Basile , laminar velocity profile. 

The theoretical analyses have generally been limited 
to infinitesimal disturbances for which the differential 
equations of the disturbance can be linearized. 

The basic problem of the small-disturbance theory is 
to determine whether a small disturbance is amplified or 
damped. If any small disturbance is amplified, it is assumed 
that a transition to turbulent flow will ultimately occur. 
However , the small disturbance theory cannot reveal any of 
the details of the actual transition because it applies 
only while the disturbance is small . The question of how 
'small' the disturbance must be in order to apply the line- 
arized equations has not been adequately answered. 

An axisymmetric flow is always stable to a disturbance 


which has only a tangential component of velocity, [Ref. 7]. 


ed 


Therefore the theoretical analyses have generally been 
limited to ASlatametr ied disturbances . gion analyses have 
been carried out by many investigators, [Refs. 2,3 and 7] , 
altough some doubt still remains concerning the details of 
the mathematics. 

Steady flow ina tube has been found to be stable to 
small, axisymmetrical disturbances by a number of investiga- 
tors,[Refs. 5,7 and 8], and may be considered as an accepted 
eonelusion. 

Several investigators have attempted to derive a ecri- 
tical Reynolds number , for a flow subjected to external 
disturbances , after which transition to turbulence occurs. 
Leite, [Ref. 9], studied the behavior of small, axisymmetri- 
eal disturbances’ in steady flow in a tube . He found that 
they were always damped for Reynolds number up to 13000 , as 
long as their amplitude was not too large. Darling,[Ref.10], 
working with pulsating flows in tubes , found that the tran- 
Sition Reynolds number’ dropped from 2510. in steady flow to 
1500 in pulsating flow. Sarpkaya,[Ref. 11], determined expe- 
rimentally the critical Reynolds number for pulsating Poi- 
seullle flow as a function of a frequency parameter and 
velocity ratio . Gilbrech and Combs, [Ref. 12], using water 
as working fluid, found a critical Reynolds number which was 


higher than 2220 , under some conditions, for steady flow . 
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They also pointed out that earlier investigators had found 
no upper limit to the upper critical Reynolds numbers. 

These experimental results confirm the conclusion that 
steady , fully developed flow in tube is stable to small di- 
sturbances , and the critical Reynolds number is a function 
of the frequency of external disturbances’ and the physical 
properties of the fluid . In spite of this confirmation , 
transition to turbulent flow is invariably observed to occur 
at sufficiently high Reynolds number and it is necessary to 
consider this fact. 

The most obvious explanation is that transition is 
caused by finite disturbances which are not damped , even 
though smaller disturbances are . This possibility is stron- 
gly supported by the observation made by Gilbrech and Combs 
(Ref. 13], that the critical Reynolds number decreases as 
the magnitude of the disturbance increases ,. 

The theoretical analysis of finite disturbances has 
been limited by mathematical difficulties , however , Stuart 
[Ref. 14], has shown that a flow may be unstable to finite 
disturbances , while it is stable to infinitesimal ones, In 
addition, Spielberg , [Ref. 15], has pointed out that stabi- 
lity to two dimensional disturbances does not necessarily 
imply stability to three dimensional disturbances. This lat- 
ter result creates a problem in comparing accurately the 


theoretical and experimental results . 
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In the present work, it was assumed that an axisymmetri- 
eal disturbance was propagated without confirming the 
non-existance of three dimensional disturbances . 

Another explanation of the origin of turbulence in pipe 
flow may be that the transition occurs in the entrance 
region of the pipe, in which the velocity distribution is of 
the boundary-layer type . A theoretical analysis by Tatsumi, 
[Ref. 2], indicates that the flow in the entrance region is 
unstable to small disturbances. He found a point of maximum 
instability near the entrance to the tube, where the maximum 
Reynolds number for stability was 19400 . However , the cri- 
tical Reynolds number’ increases rapidly in both directions, 
so that one cannot conclude that this local instability will 
lead to transition. This is because amplification required 
to produce transition remains undetermined. 

Sparrow and Lin, [Ref. 16], pointed out, in their effort 
to determine the detailed nature of the flow development in 
the entrance region, that even for a laminar flow,the velo- 
city problem for the entrance region does not yield an exact 
solution, regardless of the shape of the duct cross section. 
Their explanation for the difficulties in the analysis was 
that the essential nonlinearity of the inertia terms which 
appear in the equation of motion, preclude linear analysis. 

From a more’ sophisticated point of view , the problem 


becomes especially complicated if we attempt to reveal the 
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nature of motion in the non-linear range of boundary layer 
instability and the onset of turbulence. 

Klebanoff and others , [Ref. 17] , who attacked the pro- 
blem in the general form , pointed out that despite the suc- 
cess of the linearized theories in revealing the nature of 
the initial stages of boundary-layer instability , there re- 
mains a deep void in the undestanding of the subsequent non- 
linear behavior and the actual breakdown of the laminar 
boundary-layer. 

Miller and Fejer, [{Ref. 18] , who studied the transition 
in Blasius-type boundary-layer produced by a free stream ha- 
ving an oscillatory component of velocity , found that the 
transition Reynolds number depends only on the amplitude of 
of the oscillations , and that the dimensionless transition 
length is a function only of the reauaenesn 

Obremski and Fejer , [Ref. 19] , continued the previous 
work and pointed out that below a critical value of the pa- 
rameter [ (Re) yg LAU/2nv ] , transition occurs at a relati- 
vely constant Reynolds number , which appeared’ independent 
of the amplitude and frequency of the oscillation , at least 
over the range investigated. 

Finally Landahl , [Ref. 25], used kinematic wave theory 
to determine under what conditions breakdown of a steady or 


unsteady laminar flow into high frequency oscillations 


29 


should oeceur . Isolating the three different ingredients of 
the mechanism of breakdown , he pointed out: 


a. The primary instability process itself is only 
incidental to breakdown and transition , since one can 
easily produce a breakdown condition by a sufficient- 
ly intense local disturbance of the shear flow , for 
example through free stream disturbances. 


b. If the breakdown is to be self maintained as it 
travels downstream however , hydrodynamic instability 
is required. In a sense , the breakdown mechanism is 
the one most essential in transition of a boundary- 
layer to turbulence rather than the eclassical hydro- 
dynamie instability one, since the former represents 
a strongly irreversible process. 


ec. An unstable small amplitude wave packet of the 
Tolmien-Schlichting type may amplify for a while , as 
it moves downstream , but once it has passed through 
some streamwise position in the boundary-layer it 
will leave the shear flow practically as undisturbed 
as before the passage of the packet , except for a 
small deformation of the mean flow of second order 
in disturbance amplitude . At breakdown on the cont- 
rary there is an irreversible redistribution of the 
basie shear flow vorticity due to a nonlinear recti- 
fication mechanism . 
In view of the difficulties encountered in applying 
the small disturbance theory to steady axisymmetrical flow , 
it is likely that considerable time and effort will be re- 
quired to obtain an exact solution for unsteady flow pro- 
blems , especially if we consider a three dimensional flow . 
The use of a ' quasi-steady ' analysis would reduce the 
problem to one which ecould be handled more easily. Using 


this technique one would analyze an instantaneous’ velocity 


profile as though it were a steady-flow problem. 
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By considering the stability of the velocity profiles 
at various times, one might then draw some conclusions about 
the stability of the unsteady flow an approach used by Shen, 
[Ref. 21], as well as Greenspan and Benney, [Ref. 22]. 

In the present work, although Reynolds Numbers as high 
as 19400 were not reached , we recognize that the flow in 
the entrance region has , in all probability , affected the 
downstream flow and may well have introduced turbulence sub- 
sequently damped along the tube , for low Reynolds numbers , 


creating a laminar or laminar-like flow at the tube exit. 


B. NONSTEADY FLOW IN A CIRCULAR TUBE 

In order to study how the pressure disturbances, intro- 
duced at the entrance of a long tube used to form the jet , 
affect the flow at the exit of the tube , which constitutes 
the initial conditions of the flow of interest , we may exa- 
mine the nonsteady flow in a long tube of circular’ cross 
section . We adopt, as shown in Fig. 4 , cylindrical coordi- 
nates whose x-axis is identified with the center line of the 
pipe. 

The equation of continuity is: 


V(pV) ol = 0 
P So | (Eq. 1) 


In cylindrical coordinates Eq. 1 becomes: 


a0 


1 0 1 9 5 Ola 
0 0p 0 08 ox 





Figure 4. Cylindrical Coordinates in a Round Pipe. 





Assuming incompressibility, (o = constant), and a non- 

twisting flow on = 0), Eq. 2 may be written as: 

1 9 ; » 

—-— (rov,) +— (ov,) = 0 

eae 2 Se 

ov OV Vv 
or: eS eee ee 0 (Eq. 3) 

ox or r 


Neglecting external forces , (i.e gravity) , the Navier- 


Stokes equations in two components can be written: 








ae + V ee + Vv O¥x = aes + I Ys + ar + x _x 
ot “Ox "Or pox | - pilose Or? p UE 
(Eq. 4) 
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If we assume quasi-steady flow, the equation of conti- 


nuity Eq. 3 becomes: 


* = 0 or: Vee viret) (Eq. 6) 





This indicates that the velocity in the direction pa- 
rallel to the center line is a constant at each moment. 
Substituting an 0 into Eq. 59 we get: 


oe 0) or: P = P(x,t) (Eq. 7) 


Inserting these results in the equation of motion 








2y 1 
oYx = Se ad + |: = += —3| (Eq. 8) 
ot p x dr2 ror 
where: V= = 
Pp 


1s the kinematic viscocity of the fluid. 
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C. SOLUTION FOR A UNIFORM PULSATING FLOW 

Assuming that the pressure pulses’ introduced by the 
drive unit, (a loudspeaker driver), are propagated instantly 
along the f low direction , ( justified by the fact that the 
mean flow velocity is very small compared with the speed of 
sound in the fluid ) , with the aid of Eqs. 6, 7 and 8 we 
can show that the pressure gradient , (0P/0Ox ), becomes a 
function only of time , t. 

Using Fourier series we can express the pressure gra- 
dient as follows ;: 


coo 


o S S 
"“=— = & + =] oicosnt) +) 7 be sinnuy (Eq. 9) 
or, for the sake of simplicity of calculation , in a complex 
form: 
1 OP Ba 
. =a > int = 
> % * n=1¥n® *¥n = % 7 1D, vatele 0) 
O Ox 
where ; O°? DT are constants representing the amplitudes of 
elemental vibrations. 
The corresponding solution for the axial velocity may 


be assumed to be of the form : 


a CO 


Yx = Vx,+ = Vn, cost) + =, Vn, Sint) (Eq. 11) 
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Or, in complex form: 


Ga 
~  — = int 
[n=l ‘ne 


o> (Eqs ds) 


V = ° 


where coefficients a and mate assumed to be functions 
of the radial coordinate , r , only. 

Substituting Eqs. 10 and 12 into Eq. 8 and comparing 
the terms of the same order , we obtain the following Dif- 


ferential Equations for the coefficients as independant 


variables. 


Z 
d Vy, Ldvy, Vx 














+— +—- = 0 
dr2 r dr Vv ‘oni awt 
and: 
d2y 1 dv -in % 
en cos e+ ae repeals = 0 (Eqn 14) 
dr2 r dr v v 


After integration we get: 


ee: 
= of (E 15) 
aS sl Aj int) + By ole 
vy 4 


where: A BO = constant 


vo = - + Drcori! -) at E Ko(Srit/?) (Eq. 16) 
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where: 98 = V(n/v) 
and Jo and Ko are the Bessel functions of the first and 
second kind of zero order respectively. 

The solution for the x-wise component of velocity is 


given by the real part of the following expression: 


% 2 
x Vv 4 O O 


oS 
+ =, =. i Dado (Sri?!) fe E Ko(Srit/?) eint (Eq. 17) 


On the axis of the pipe , ( r=0 ) , In(r) = -° and 
Ko(Srit/*) = 00 » Since vy, is finite , we conclude that: 
i E = 0 (Eq. 18) 


The Boundary Condition at the wall is given Dy: 


- 


v. = 0 at re=R (Eq. 19) 
wheres R is the radius of the pipe. 

Inserting this boundary condition , Eq. 19 , into 
Eq. 17 we get: 


Ds oo 


- aR Ae 
5 een= 1 


2 4v in 


3/2) 


+ DJ, (8Ri (Eq. 20) 





Inserting , Eqs. 18 and 20 into Eq. 17, we finally get 


the velocity component Vy ° 
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ain ne 
Vv see [R2-r2] = on 1 - 2 aint \Ea: 21) 
n J (ORI ) 


The total mean mass flow , m , which is transmitted in 


the x-direction is given by the integral over time and ra- 


dius of the axial velocity component: 


—_— CEGe 227) 
é } 8v 


Inserting the time mean pressure gradient: 


= (Bq. 23) 
-_ — = re4 Gene 
ox pmo 
we get: 
aR*t dP 
m= ——|-—— (Eq. 24) 
8y Ox 
lt 


is seen that the total mean mass flow , in pulsating 


motion , is identical to that of a steady Poisseuille flow 


with the same pressure gradient, 


as the mean pressure gra- 
dient in the pulsating flow. 


A mean velocity , U , may be defined by: 
m 
Uz —— 


9 (Ea. 25) 
TR 
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Combining Eqs. 24 and 25 we get: 


R2 





‘OP 
0 | (Eq. 26) 
Sy Ox 


We may express the x-wise velocity as the sum of a 


e © ‘ ° 
perturbation velocity, v, . and the. mean stream velocity of 





the fluid,v,, ° 
Ve =F Vy5 * ve (Eq. Je) 
Introducing the Hagen-Poiseuille relation for Vig we 
t 
ge : 9 
XS 
U R 


Dividing Eq. 27 by the total mean velocity , U , we 


get the non-dimensional expression: 





(Eq. 293) 


Combining Eqs. 21 and 29 we get for the perturbation 


velocity , Ve : 
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where: 


ber(OR)ber(dor) + bei(OR)bei (dr) 





ber“(8R) + bei*(OR) 
bei(dR)ber(dor) - ber(dR)bei (or) 


ber“(8R) + bei2(5R) 


(Eq... 34) 





oo) 
i" 


wheres: ber and bei are the Kelvin functions . 


With the aid of Eqs. 10 and 26 we can obtain the cor- 


responding pressure gradient relation: 


oo 


= Pn sinnt)]} (Eq. 32) 
n=l 
Xx 
O O 


ao 
4R OP 64 x 
= —(1 + = =! eos(nt) + 
n=1 Fy 


Introducing the friction factor for laminar flow : 
f= 64 / Re 


into Eq. 32, we get : 


oo 
a b 
<— —cosnt) + a, siapy (Eq. 33) 
Ox 4R Oo, ays 


where: % /q, and bn/do are dimensionless ratios of 


amplitude of periodic variation of pressure gradient norma- 


lized with the mean amplitude . 


395 


D. VELOCITY DISTRIBUTION AND ASYMPTOTIC EXPRESSIONS 

The velocity profiles for laminar , pulsating flow cons- 
titute the initial conditions for any theoretical analysis 
of stability. Although no exact theoretical analysis of sta- 
bility is presnted in this work, solution for the velocity 
profiles will serve to identify the dimensionless parameters 
which are significant and will permit some qualitative de- 
duction concerning Oye ancien ety of pulsating flow, 

From Eq. 16 we found that the periodic part of the velo- 
city profile is characterized by the parameter 5R =\/n/v R. 

The distribution of a,/ a, , the dimensionless velocity 
amplitude, obtained by Gilbrech and Combs , [Ref. 16] , is 
plotted in Figure 5, 

Uchida, [Ref.23], obtained solutions for limiting values 
of the parameter oR by means of asymptotic expansions of the 


va 


Bessel function: 


1. Slow Oscillation ( 8R << 1 ) 
When a highly viscous fluid oscillates slowly in a 


narrow pipe , OR becomes small and we can approximate: 


ber(OR) — 1 and bei(dR) — 0 


Using this approximation, we get from Eqs. 28 and 29 


2 — a — }b 

V r > Nia 

eS 2[i-—]|: + =, —ecosint) + n=l 7s ini) 
R2 - OX Xo 
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Figure 5. Oscillation Amplitude of the Local Velocity 
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2 s) % 
a r a Pe 
ors: —— = 2 =| | > | 
U Rt 
1 ' 108 
: ee 2_. r ) = =e (Eq. 34) 
or Vv, = Pe | Se 


From this we econelude that the velocity distribution 
is a parabola, as in the case of steady Hagen - Poisseuille 
flow, while the magnitude of velocity varies periodically in 


phase with that of the pressure gradient. 


2. Rapid Oscillation ( 6R —»© ) 
When a fluid of low visecocity oscillates rapidly in 
a large pipe , then parameter oR becomes large. 
Assuming 68R > 10 , asymptotie expansion of Bessel 


functions may be introduced: 


gu «4 loon) oe 
ber (oR) ask cos — . = 


eOR/2 OR 1 
bei (OR) fs Sinl|— - —| + 0 =| 
2 75R 2 8 


Eq. 35) 


To examine the motion of the fluid near the center 
of the pipe , we put: 
OR —=co and or —= 0 


As before we get: 








V r a > 
“so af-h) + Se sony - B28 —geomny 
U R2 a, (SR) a (SR) 








V 2 = o4 
r 8 TI 
ee xii 2f1- : = eos (nt -—| + 
a B a, (SR)? 
GO 
b 8 ae 
Ls = sin|nt -—| 
%, (OR) : 


ors Vv = —2 (R2_ r*) ‘ao cos|n =] + +2 —tsin[ar- =] 
4v 2 
(Eq. 36) 

From this result we may conclude that in the case of 
rapid oscillation , fluid flows at the center of pipe with 
the phase lag of 90° behind the wave of pressure gradient 
and its amplitude diminishes with increasing frequency. 

To examine motion of the fluid near the wall of the 
pipe , we puts: 

SR —= co and 6R -—co 

Using the same procedure we get: 


9 CO d(R-r) Ps) 
=== 9 ee 1 —; es sin|nt-—(R-r) 


§(R-r) 


R § 
2 jie | =. Cos pecan cea, 37) 
i co s(n t)+ mC 2 | E 


ete oF 
na (3K)? 





Because of the obvious complexity of the expression, 
it is difficult to draw any simple conclusions other than to 
note that a local maximum in velocity developes in the neigh 
borhood of the wall as it its shown in Fig. 5. We can also 
get a better confirmation from Figs. 6 , 7 and 8 , obtained 


from (Ref. 23], for values of dR = 1, 3 and 10 respectively. 
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Figure 6. Velocity Profile for 6R = 1 
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IR=3 





Figure 7. Velocity Profile for 


dR 


i 
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-0-10 


Figure 8. Velocity profile for 6R = 10 
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E. SECTIONAL MEAN VELOCITY 


Since the instantaneous mass flow and accordingly sec- 
tional mean velocity , are periodic functions » we seek to 
link the mean velocity with the applied pressure gradient. 


Sectional mean velocity , denoted by Vin, 18 given by? 


1 & 
V = ——_— y) 
om = wi Rive cr (Eq. 83) 


Introducing Eq. 34 and integrating we get: 








V sa oO 8 
_ 1 + est oan 
10 4 





2D 2C 
—cos(nt) + 1- —— sin 25 
dR 








U ee (SR)*15R 
S “s 8 2D . | 2G | 
1 ae ri ae f=) | l= Se] Cosgm 
mt ao, (SR)? = sel oo] (Ea 39) 
where: 
d[ber(dR) ] ) d[ber (dR) ) 
ber (dR) a athe , bei(OR) = —————__ (Eq. 40) 
d(oR) d(dR) 
ands 
ber(oR)bei(OR) - bei(dR)ber(dR) 
Gea 
ber“(5R) + bei“(SR) 
; (Eq.41) 
ber(OR)ber(OR) + bei(d5R)bei (dR) 
RS SSS ——————EEee 


ber“(SR) + bei’(dR) 


An example of simple periodic pulsation is given by: 
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= sce 
neosnt 


.—_ 

——y 
i 

4 

rg 

EE — 
i 


Inserting the above relation into Eq. 











( average term ) 
term ) 


( perturbation 


39 we get: 


V om g [2D AC 
iw a —cosnt) + F ‘ | sine 
U a (OR), LOR OR 
Vin a, 8 2c}? 2D)? 
ONG a 2 ie Se 5 cerca (a: | cos(nt ~ $,) 
U o, (8R) 5R OR 
v a 
xm n 
Or: ———= 1 +-—A_ eos(nt - >) 
U a * " (Eq. 43) 
Oo 
wheres 
A = 
V 





is the coefficient of amplitude , and: 


1-2C/65R 





2D/0R 





OR-2C 








2D 
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represents the coefficient of phase lag in the wave of 
pressure gradient. 

In Fig. 9 , obtained from [Ref. 23] , are plotted the 
coefficients of amplitude and phase lag of mean velocity as 


a function of the parameter oR, 





Figure 9. Coefficients of Amplitude and Phase 


Lag of Mean Velocity. 


F. SURFACE FRICTION AND BALANCE OF FORCE 


The instantaneous frictional force acting on the pipe 


wall is given by: 


— 
t= -pl|— 
P dr 


(Eq. 44) 





ar | 
i 
a 


Dividing Eq. 44 by ( pU*/2 ) , and combining with the 


previously obtained Eqs. 32 , 39 and 40 we can get: 
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Zt 1 2D 
——_ = mae a a | — —cos(nt) + —sinnd| i 
a OR 


pu? Re 5 OR 
— ey UC 2D 
eee ee t (Eq. 45) 
n=l — | oe cosh 


O 
Considering the simple periodic pulsation described by 


Eq. 42 we get: 

















7 16 a, IG 2D 

—— = ——/1 + | —cos(nt) + —sint) 

Hag 16 x, 
or: = a eC O Sian ee 

9 4E 
pU Re re 
4ovU 

or: vas 1 neal qeeos(nt - >.) (Eq. 46) 
where: “o 





2 2 


2D 
OR 


2C 
SR 














represents the coefficient of amplitude of shearing stress 


at the wall , and: 


?_ = tant] —] 
C 


represents the ecoefficient of phase lag of shearing stress 


with respect to the wave of pressure gradient. 
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In the Fig. 10 below , obtained from [Ref. 23], it is 


shown how the coefficient of amplitude and phase lag of 


shear stress varies with 5R 





Figure 10. Coefficients of Amplitude and Phase Lag 
of the Shearing Stress 

If we compare Figs. 9 and 10 we observe no substantial 
differences in amplitude or phase lags. Moreover, the phase 
of the sectional mean velocity is greatly delayed with res- 
pect to the pulsating wave of pressure gradient, while that 
of shearing stress is less delayed from it. 

Using Eq. 8 and integrating over the section of the 
pipe, from r=0 to r=R , we can obtain a relationship between 
the force and acceleration of amass of fluid enclosed in the 


circular cylinder of unit length: 





OP dv 
- m1R2— = R“p + 2 RT 
ox dt 
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LG oP/ox ; dv_,/dt 4 


or: (Eq. 47) 








ees i 
2 pu" /A4R U*/2R pu" 


For the present example of a simple periodic pulsation 
we find the mean frictionel force , involved in the shear 
stress , is just balanced by the term of mean pressure gra- 
dient and there results no steady acceleration. The various 
terms appearing in Eq. 47 are shown in Fig. 11 below , ob- 


tained from [Ref. 23]. 


Az Re -dP/dx 
ABC 32 pV*/2R 


Re 2% 
32 pv¥2 
>  & iis 
gp. Fe tum /dt 
32. V*/2R 


Figure 11, Periodic Force Component (A - Pressure Gradient, 
B - Acceleration , C - Shearing Force) 
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IV. RESULTS AND CONCLUSIONS 


A. PRELIMINARY RESULTS 

After a preliminary survey of the phenomenon of jet 
breakdown at a single laminar and single turbulent Reynolds 
number, as mentioned above, it was found that the field was 
so rich in information that a third , transitional Reynolds 
number of 2200 was added to the survey . Observations repor- 
ted are of the most dinstinetive phenomena at each of the 
three Reynolds numbers and were observed to occur most fre- 
quentlyat frequencies which were multipliles of 33 and 50 hz 
Results inelude two photographs at each data point of two 
separate phases along with an oscillogram showing both the 
excitation waveform and the pressure perturbation it produ- 
ced, shown in Figs. 13 through 36. Table I is a summary of 
the observed breakdown length as a function of the frequency 
and Table II is a summary of oscilloscope and flow parame- 
ters, corresponding to each figure. 

As Hoyt and Taylor, [Ref. 24], pointed out, jets are un- 
stable to two kinds of disturbances, axial and helical. The- 
ory suggests that the axial instability ( corresponding to 
the initial waves on the jet surface ) die out after a few 
nozzle diameters downstream and the helical instabilities, 


amplified by aerodynamic resistance , lead to the final jet 
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breakup. The initial slight helical motion is amplified un- 
til the entire jet may be described as having a ecorksrew- 
like motion . This is well confirmed by the photos obtained 
in the present work of the jet formation before and after 
the breakdown and is in good agreement with the observations 


obtained by Crow and Champagne, [Ref. 4]. 


TABLE | 


SUMMARY OF BREAKDOWN DISTANCE 


f (Hz) #d #d #d 
(Re=1760) (Re=2200) (Re=2640) 

0 118.0 125.0 30.0 
33 112.6 120.6 31.9 
20 110.0 112.0 34.0 
66 106.3 103.7 32.2 
100 102.0 98.5 29.9 
133 98.6 99.2 33.7 
150 95.95 100.0 35.0 
166 92.7 101.7 37.8 
200 88.0 102.5 41.0 
233 85.2 101.2 49.2 
250 83.5 939.0 49.0 
266 81.95 101.5 47.6 
300 91.2 88.0 45.0 
333 96.1 97.5 44.1 
350 101.8 110.3 43.3 
366 110.0 121.0 42.0 
400 112.0 125.0 42.9 


Re = Reynolds Number 
f = Introduced Frequency 
#d = Number of Diameters 
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TABLE II 


SUMMARY OF EXPERIMENTAL CONDITIONS 


Fig.# Re f V 
p Vg t /em 
(Hz) (mV/em) (mV/em) (msec/em) 


13 1760 90 =) 10 a) 
14 1760 100 ") 10 -) 
15 1760 150 ") 10 2 
16 1760 200 ) 10 2 
17 1760 290 ) 10 2 
18 1760 270 =) 10 2 
19 1760 300 =) 10 2 
20 1760 380 ") 10 2 
21 2200 20 =) 10 ") 
22 2200 100 ") 10 ) 
23 2200 150 ") 10 2 
24 2200 200 =) 10 2 
29 2200 233 =) 10 2 
26 2200 250 ") 10 2 
27 2200 266 ) 10 2 
28 2200 300 =) 10 2 
29 2200 380 10 10 1 
30 2200 787 20 20 0.9 : 
31 2640 30 ") 10 =) 
32 2640 100 =) 10 2 
33 2640 200 *) 10 2 
34 2640 250 10 10 2 
35 2640 300 10 10 1 
36 2640 600 10 10 0.9 

Re : Reynolds number 

f : Introduced Frequency 

Mis Signal Amplitude from Pressure Transducer 

V 


g : Signal Amplitude from Frequency Generator 
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Figure 12. Breakdown Length Versus Frequency 
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1. Laminar Flow at Reynolds Number 1760 


At this Reynolds number’ the unexcited flow started 
to breakdown at approximately 118 diameters from the tube 
exit. Increasing the frequency to 266 Hz produced a gradual- 
ly movement of the breakdown point upstream to approximately 
81.5 diameters from the exit. Above this frequency a rever- 
sed trend appeared and the breakdown point started to move 
downstream as the flow became laminar similar to the unexci- 
ted flow , the only difference being in the droplets formed 
after the breakdown. This phenomenon may be characterized as 
a relaminarization of the flow and was almost complete at 
380 Hz. Further increasing of the frequency up to 1000 Hz 
produced no significant effect in the jet flow. 

Results for Reynolds number’ 1760 are shown in Figs. 
13 through 20 . As can be observed from the photos , break- 
down point is not fixed and oscillates somewhat upstream and 
downstream , creating a breakdown region . The amplitude of 
this oscillation was reduced as the perturbation frequency 
was increased , starting from about 14 diameters and approa- 
approching 4 diameters as the frequency changed from 200 to 
266 Hz . At 266 Hz an especially interesting phenomenon was 
observed concerning the breakdown point and the droplet for- 
mation . The breakdown of ane jet occured at the minimum ob- 


served distance from the tube exit at 81.5 diameters. 
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Photos Illustrating Jet Breakdown and Oscilloscope 


Figure 14. 
Trace of Pressure at Re=1760 and f=100 Hz 
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Figure 17. Photos Illustrating Jet Breakdown and Oscilloscope 
Trace of Pressure at Re=1760 and f=250 Hz 
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Figure 18. Photos Illustrating Jet Breakdown and Oscilloscope 
Trace of Pressure at Re=1760 and f=266 Hz 
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Figure 19. Photos Illustrating Jet Breakdown and Oscilloscope 
Trace of Pressure at Re=1760 and f£=300 Hz 
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Figure 20. Photos Illustrating Jet Breakdown and Oscilloscope 
Trace of Pressure at Re=1760 and f=380 Hz 
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The oscillation of the breakdown point was minimized to a- 
bout 4 diameters and the droplets formed , as was observed 
with the aid of the stroboscope , were exactly the same and 
appeared stationary in space. This is shown in the jet pho- 
tos of Fig. 18. 

A plot of the breakdown distance as a function of 
‘frequency is provided in Fig. 13 where we see that, in the 
frequency range from zero to 266 Hz, an almost linear varia- 
tion of the breakdown point with frequency occurs. The mini- 
mum at 266 Hz may be characterized as a eritiecal frequency 
for this speecifie Reynolds number . Above this frequency one 
also observes an almost linear variation of the average 
breakdown point to 380 Hz. This region may be characterized 
as relaminarization region. The slope of the breakdown curve 
is greater than at lower frequencies although increasing the 
frequency above 380 Hz results a zero slope curve , as 
discussed above. 


2. Transitional Flow at Reynolds Number 2200 


At this Reynolds number the unexcited flow started 
to breakdown at approximately 125 diameters from the tube 
exit . Increasing the frequency to 100 Hz , the breakdown 
point moved upstream up to about 98 diameters from the 
tube exit . The oscillation of the breakdown point was 


reduced at this case, from 12 to 6 diameters. 
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Increasing the frequency to 266 Hz resulted in 
only slight variation of the breakdown point . The only 
unusual observation in this range of frequency was the 
drastic reduction of the amplitude of breakdown point 
oscillation as the frequency was increased ,. 

An interesting pressure effect was observed at 266 
HZ, as can be seen at Fig. 27 , where a maximum amplitude 
of pressure signal was recorded. 

An interesting change was observed when the fre- 
quency reached 300 Hz . The breakdown point was displaced 
suddenly about 25 diameters upstream and nearly stabilized 
at 88 diameters . The droplets viewed with the aid of the 
stroboscope appeared immobile and any oscillation of the 
breakdown point was difficult to observe . Based on a series 
of high speed photos it was concluded that the breakdown 
point was oscillating with an amplitude of about 2 diameters 
Also observed at this frequency that the waveforms obtained 
from the Signal generator and the pressure transducer’ were 
almost in phase and moreover, the latter waveform appeared 
nearly sinusoidal as shown in Fig. 28. 

Further frequency increases resulted in a flow rela- 
minarization as was observed and discussed in the previous 
case. A completely laminar flow was observed at a frequency 
of 380 Hz and no additional significant effects were obser- 


ved up to 2000 Hz . 
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Trace of Pressure at Re 


Figure 21. 
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Photos Illustrating Jet Breakdown and Oscilloscope 


Pigwre 22. 
; Trace of Pressure at Re=2200 and f=100 Hz 
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Figure 23. Photos Illustrating Jet Bre 
Trace of Pressure at Re=2200 and f=150 Hz 
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Figure 24. Photos Illustrating Jet Breakdown and Oscilloscope 
Trace of Pressure at Re=2200 and f=200 Hz 
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Figure 25, 
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Figure 27. Photos Illustrating Jet Breakdown and Oscilloscope 
Trace of Pressure at Re=2200 and f=266 Hz 
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Figure 30. Photos Illustrating Jet Breakdown and Oscilloscope 
Trace of Pressure at Re=2200 and f=787 Hz 
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Refering to Fig. 12 , we may draw some conclusions 
eoncerning the average breakdown point. Up to 100 Hz a gra- 
dually movement of the breakdown point upstream is observed 
which is nevertheless more rapid than the laminar case as 
ean be seen from comparison of the_ slope of the curves | 
within this region . A nearly stable region then appears in 
the frequency range between 100 Hz and 266 Hz where the 
breakdown point is nearly fixed at 100 diameters. This re- 
gion may be characterized as a_ stabilized region and was 
not observed at any other Reynolds number. A_ singular point 
was also observed at 300 Hz as discussed previously. This 
frequency may be characterized as a critical frequency for 
the specific Reynolds number . Finally , we can see _ from 
Fig. 12 that the relaminarization rate was also greater for 
transitional flow, Mrecuc tie a steeper slope of the curve in 


the frequency range from 300 Hz to 400 Hz. 
3. Turbulent Flow at Reynolds Number 2640 


At this Reynolds number it was’ observed that the 
unexcited flow began to breakdown at an average distance of 
30 diameters and no significant effect on the jet was obser- 
rved up to a frequency of 100 Hz. Increasing the frequency a 
gradual movement of the average breakdown point downstream 


observed, which reached a maximum of 50 diameters at 250 Hz. 
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Further increase of the frequency produced a very 


slow movement of the average breakdown point upstream to 


about 42 diameters at a frequency of 366 Hz. The breakdown 


point then remained almost unchanged up to 2000 Hz . The 


only further change observed in the high frequency range , 


from 


366 Hz to 2000 Hz, was the appearance of a corksrew- 


like motion of the jet which became more’ pronounced with 


increasing frequency , making more readily observable the 


helical mode , which dominates the breakdown mechanism , 


These 


results are shown in Figs. 31 through 36. 


B. CONCLUSIONS. 


From the results , discussed above , we may draw the 


following conclusions concerning the effect of small pertur- 


bations on jet breakdown: 


a. 


In the case of laminar flow, (Re=1640) , the breakdown 
point moves upstream as the frequency is increased up 
to 266 Hz. Downstream of the breakdown point the jet 
flow becomes’ turbulent and breaks’ into discrete dro- 
plets. Above 266 Hz this effect is observed to reverse 
itself and the jet returns to its laminar nature. 


In the case of transitional flow , (Re=2200) , the 
breakdown point moves initially upstream as frequency 
is increased to 100 Hz . In the range of 100 Hz to 
266 HZ an almost stable behavior of the breakdown 
point is observed . Singular behavior with very inte- 
resting results including a sudden shift of breakdown 

point is observed at 300 Hz . Finally a more effective 
relaminarization than observed at lower Reynolds num- 

ber occurs above 300 Hz leading to a completely stabi- 
lized laminar flow at frequencies above 380 Hz . 


In the case of turbulent flow , (Re=2640) , no notice- 
able effect is observed up to a frequency of 100 Hz. 
Above this frequency the breakdown point begins to 

move gradually downstream up to 250 Hz , returns back 
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upstream to a frequency of 333 Hz , and becomes almost 
completely unaffected by the pressure perturbations 
at higher frequencies. 


Changing the waveform of the disturbance produced 
only a slight deviation inthe average breakdown 
point, however a markedly different patern of droplet 
formation was observed. 


Although it was not possible to record the phenomena 
photographically for inelusion here , it was possible 
using the phase angle control incorporated in the 
stroboscope to observe the entire development of the 
droplet formation visually. This process was found to 
be a complex sequence of breakup and re-agglomeration. 
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Figure 31. Photos Illustrating Jet Breakdown and Oscilloscope 
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Figure 32 


RAYA 


- 

— 

he 
> 





rr tae = 
3 4 
? @ 
o 
ro) 
[4 
2 2 
ca ee 
> 4 
° 
a 
~ at2tihSoeue 
= 
€e 
: oe 
r 2 , > 
rs : 
: ; 4 
eo = 
q a z a 
arg S 
“ - 
, ‘ ‘ 
een : 
° 
¢ one 
° 
4 5 
* 
* 
4 ; 
ria = 
, e 
peers cowie: = 
eal 
: . 
~ | oe j 
° 
—_ 
abe 
a 
* 


oat ees aS Sees cee 


Figure 33. Photos Illustrating Jet Breakdown and Oscilloscope 
Trace of Pressure at Re=2640 and f=200 Hz 
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Figure 34, Photos Illustrating Jet Breakdown and Oscilloscope 
Trace of Pressure at Re=2640 and f=250 Hz 
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Figure 35. Photos Illustrating Jet Breakdown and Oscilloscope 
Trace of Pressure at Re=2640 and f=300 Hz 
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_ APPENDIX. CALCULATIONS AND DATA REDUCTION 


A. CALCULATION OF CRITICAL REYNOLDS NUMBER 

A visual observation of the flow was used in order to 
determine the critical See noid number . Adjusting the flow 
Dyethe means of the needle valve , installed at the entrance 
of the plenum chamber , a laminar jet was observed up to the 
flow meter indication m, = 25 . Above this value a turbulent 
flow resulted. 

Using the equation of continuity 

m= pAU or U=m/ pA (EaiesA— 1 ) 
where: 

m = the fluid mass rate. 

U = the mean flow velocity. 

9 = 62.34 lbm /ft® , the water density (at 60° F). 

A = the cross section area of the pipe. 


a Gree 16) 129) (et?) 


or: Av= sieaa7s yao? (ft?) 


A = trd 
and: 
m = [.6(25/100) J (gal/min) or m = .15 gal/min 


or: m = (.15)(500.8)(1/3600) (lbm/sec) 


or: m = .020867 (lbm/sec) 


Substituting the above values into Eq. A-l we get: 
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. 020867 
Ue (ft/sec) 
(62.34)(1.9175 1074 





or: U = 1.74562 (ft/sec) 
The Reynolds number is given by: 
Re = (pUd) / (pg,) GEG A~2) 
y = 2.4 .107° (Ibf-see/ft*) 
and: 
g, = 32.2 (1bm-ft)/(Ibf-see*) 


Substituting the above values into (Eq. A-2) we get: 


(62.34)(1.74562)(1/64) 


(Re) 
a (2.4 107°) (32.2) 
Or: (Re). = 2200 | 
The value obtained may be considered reasonable taking 
into account that is the most common value found in the Hy- 
drodynamics Bibliography. The numerical values’ used in 


the above calculations were obtained from Ref. 25 , 


B. CALCULATION OF INDICATED PRESSURE IN LAMINAR FLOW 
Refering to Fig. B-1, where the exact dimensions used 


are indicated, we apply the Energy Equation: 





g 
+2, = ae zo + hy (Eq. B-1) 
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DETAIL A 





| ; 
24" 
| 36° d= 0.1875 
¥ —— | 
2,715) d,=0.25 
Se . —— REFERENCE PLANE 
DETAIL A 
Figure B-1. Sketch of the Apparatus. 
where: 
P, : the pressure at point l. 
P, ; the pressure at point 2 ( P, = Bein = 14.7 psia) 
U, ; the mean velocity of the fluid at point l. 
U, the mean velocity of the fluid at point 2. 
bs = be = 2 : the energy correction factor for laminar 


flow at points 1 and 2 respectively. 





Zi1Z9 3 the height difference from the reference plane 
ny : the flow losses given by the relation: 
aU" 
hy = i == (Eq. B-2) 
d 28 
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where: 
f = 64/Re : friction factor for laminar flow. 
I :; the total pipe length. 
Inserting the numerical values and the appropriate di- 
mensions from Fig. B-1 into Eq. B-2 we get: | 


(64/2200)(94.416/.1875) (1.74562) “/2(32.2)] 


> 
i 


ors h, = .69335 (lbf-ft/lbm) 
Using the Continuity Equation refering to Detail A of 
Fig. B-1l and neglecting the minor turning losses we get: 
m = pAU, = pAU, (Eq. B-3) 
: = zs = 2 aa 
or: U,/U, = A,/A, = d do = (.25)°/(.18%5)> =" 77 


Therefore: 


Cc 
tI 


U,/1.7777 = 1.74562/1.7777 Cfit/see) 


982 (ft/sec) 


© 
ar | 
on 
— 
t 


From the Energy Equation, (Eq. B-1) we may write: 


4 b 
—~ (Py = Pe) ~ 2 (U,?-U,7) + Blzp-24) + Behy 


p 
Solving for AP = Pi-P, we get; 
PL (0,202) aaa! 
AP = — (— (U,4-U + g(Z,-z,) + gah 
g. 2 2 1 2 %1 con (Eq. B-4) 


Substituting numerical values into Eq. B-4 we get: 


62.34 
Ape a ee (1.74562) 2-(.982)2]+32.248+32.2(.69335)} 
32.20% * le 


or: AP = 125.18 (lbf/ft*) 


86 


or: AP .8693 (psig) 


or: AP = 1.775 (in.Hg) at 60° F 


The latter is a reasonable result compared with the 
mercury manometer which indicated pressure P sad = 1.8 in.Hg 
suggesting an error of 1.39 % . 

The curvature of the pipe does not create a large effect 


in friction faetor for this Reynolds number’ which is in 


good agreement with the experimental work of White [Ref.26]. 


C. DYNAMIC RESPONSE OF INDICATED PRESSURE 
Since the pressure perturbations were introduced in 
the upstream plenum chamber and measured with a pressure 
transducer through a pressure tap, we may reasonably inquire 
into the accuracy of the pressure traces measured. In order 
to answer this question, we consider the dynamic response of 
the system under the assumption that the fluid pressure is 
applied uniformly at the transducer front surface . 
The transient response of a pressure measuring instru- 
ment is dependent on two factors: 
a. The response of the transducer element that senses 
ses the applied pressure. 
b. The response of the pressure transmitting fluid 


contained in the connecting tubing. 
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This latter factor is frequently the one which deter- 
mines the overall frequency response of a pressure measure- 
ment system , and eventually calibration must be relied upon 
for determining this response. 

An estimate of this behavior may be obtained anor the 
following analysis [Ref. 27]. 

Consider the system shown in Fig. C-1. The fluctuating 


pressure has a frequency, , and an amplitude, P and is 


0’ 
impressed on a tube of length L, and radius r. At the end of 
this tube is a chamber of volume V, where connection with 


the pressure transducer is made . We seek an expression for 


the pressure indicated by the transducer, Poe 







P, sinut 


| PRESSURE TRANSDUCER 
—_ | 


Figure C-1. Sehematiec Diagram of Pressure 


Transducer Connection. 


The mass of the fluid oscillates under the influence 
of fluid friction in the tube which tends to dampen the mo- 
tion. If the conventional formula for laminar frietion resi- 
stance in tube flow is used to represent this friction , the 


resulting expression for the pressure amplitude ratio is: 
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lage = {[1 - |" + an? (open? 1/2 (Eq. C-1) 
0 
where: 


Wo is the natural frequency given by the expression: 


rc”. 


oe a (HaqcnC-2) 
n ALV 


and h is the damping ratio given by the expression: 


Ee Shab | Ufa), (Cas) 
per? \ tr 


In the above relations, Eqs. C-2 and C-3 





ec = the velocity of sound in the fluid. 


the dynamic viscosity of the fluid. 


5= 
it 


p = the fluid density. 
The phase angle for the pressure signal is given by: 


Ww 
ee eee (Eq. C-4) 


Ww )2 
Le erie 
Introducing the physical dimensions of the present appa- 


ratus, refering to Fig. C-1 we have: 


Peace o. SD 1n tee Oso 1n 
fo Ey 8 VAS el r’= 0.25 in 
7 ae y) _ -5 3 
Nee=ennmeiy 8411 425/12) (3.5/2) = 6.8177 10 (ft*) 
en= (K/o)+/? 
is 3 Wo 
Orsmemce=mt44 (311 .10° )/1.938] = 4800 (ft/sec) 
9 = 62.34 lbm/ft® 
p= 2.4 .10°° lbf-sec/ft® (at 60°F 
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Substituting the above values into Eq. C-2 we get: 


37r(.25/12)2(4800) 2 


_ | _ 
Wn 4(3.5/12)(6.8177 1074) 
or: ww, = 10885 secr- 


Substituting into Eq. C-3 we get: 


2(2.4 1079)(32.2) 






3(3.5/12)(6.8177 1074) 


SNS ees 








(62.34)(4800)(.25/12) 1 


7.87 10° 


On: h 
From these results we ee that the difference between 

pressure transducer indication and actual plenum chamber 

pressure is negligible for the working frequency range. 


For example at W = 1000 Hz , from Eq. C-1 we get: 


P 
——— l- 
P 
O 


Ore P/P, = 1.0085 


1000 100012\-1/2 








212 
| + 4(7.87 1078)2 











10885 10885 


Therefore the pressure indication error due to dynamic 
response is 0.85 % 


We can get for the phase angle of the pressure signal, 


from Eq. C-4 : 


- 
tl 


tan + -[2(7.87 10-°)(.09187) / (1 - .09187)] 


=3.355 LOmeerad 


or: & 
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D. CALIBRATION DIAGRAMS 


In order to compare theoretical and experimental re- 
sults concerning the indicated pressure from the mercury ma- 
nometer, as well as the pressure variations from the digital 
multimeter, it was necessary to calibrate both of them. 

A pressure calibration diagram is plotted in Fig. D-l where 
it is shown to be linear. 

Since the Reynolds number is a function of flow rate, it 
was easier to adjust the flow more accurately with the aid 
of a pressure indicating digital multimeter than the rotame- 
ter. Therefore two diagrams were created, shown in Figs.D-2 
and D-3 where are plotted the indicated flow rate versus the 
indicated pressure , expressed in inches Hg and D.C Volts 


respeerively. 


on 
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Figure D-1. Pressure Calibration Diagram 
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FLOW RATE VS P(IN.HG) 
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Figure D-2. Flow Rate vs. Pressure Indicated (in.Hg) Diagram 
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FLOW RATE VS P(D.C VOLT) 


60.0 70.0 60.0 


40.0 


FLOW RATE INDICATED 
60.0 





6.8 7. 


6.4 6.6 0 
PRESSURE INDICATED (D.C VOLT) 


Figure D-3. Flow rate vs. Pressure Indicated (V.DC) Diagram 
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